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Abstract 



In this paper we argue that the pattern of cell movements in the morpho- 
genetic process known as convergent extension can be understood as an energy 
minimization process, provided the cell-cell adhesive energy has a certain type 
of anisotropy. This single simple property suffices to cause the cell elongation, 
alignment, and intercalation of a cellular array that are the characteristics of 
convergent extension. We describe the type of anisotropy required. We show 
that the final aspect ratio of the array of cells is independent of the initial 
configuration and of the degree of cell elongation. We find how it depends on 
the anisotropy. 
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FIG. 1. Intercalation. Isodiametric cells (a) become elongated and aligned (b)while simulta- 
neously intercalating (c) so that an array of cells extends at right angles to the direction of cell 
motion. 

In the development of the animal embryo great changes of form (morphogenesis) take 
place jlj], especially during gastrulation when axial structures are formed by extensive cell 
rearrangement. During these rearrangements groups of cells move coherently over distances 
very large compared to cell dimensions. This process has been extensively investigated by 
experiments, particularly on embryos of the frog, Xenopus laevis, and particularly by R.E. 
Keller and his collaborators (see || for a brief review and extensive references). 

One characteristic and widespread type of rearrangement, "convergent extension", oc- 
curs, for example, in the development of axial structure such as precursors to the vertebrate 
spinal column. Here an active group of cells undergoes a threefold process. The individual 
cells, originally roughly isodiametric (Fig. |I]a), elongate and their axes of elongation be- 
come aligned. If these were the only motions the final configuration would be as in Fig. 
|l]b. But at the same time, though on a somewhat slower time scale, the cells intercalate 
between each other. The intercalation is in the direction of alignment so that the number of 
cells in that direction decreases while the number of cells in directions perpendicular to the 
alignment increases producing a final configuration as in Fig. |I|c. The elongation increases 
the overall length of the group of cells in the direction of alignment and decreases the length 
in orthogonal directions (since the volume stays roughly constant). Intercalation does the 
reverse but dominates, so that the axis of net extension of the group of cells is at right 
angles to the axis of individual cell elongation. Here we argue that the important aspects 
of convergent extension will result from a tendency of the active cells to minimize their total 
energy, provided that they interact with a non-uniform surface (adhesive) energy satisfying 
certain conditions. We also develop a mean field theory of this process. 

Minimum energy principles have been used to explain cell rearrangement since Steinberg's 
H suggestion that differential cell adhesion plus cell motility can account for cell sorting 
patterns in mixtures of two or more cell types (see M for a review and extensive references to 
the literature). Goel and Lieth || have considered cell sorting for a simple geometrical model 
in the presence of anisotropic surface adhesion between cells of fixed shape. Cell sorting, 
driven by energy minimization, has also been the subject of many computer simulations 
|||7|]. Drasdo, Kree, and McCaskill M have simulated cell sorting with anisotropic surface 
adhesion. Anisotropic surface adhesion has not, so far, been used to explain the convergent 
extension of a homogeneous group of cells. We do not model here the dynamics of convergent 
extension. We assume, as in and j5|, that cell motility will allow the system to explore 
its possible configurations and that, as a strongly dissipative system, it will evolve towards 
the configuration of minimum energy. 

In the embryo convergent extension usually takes place in an asymmetric environment 
where the inactive cells bounding the active region are not the same on all sides of that 
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FIG. 2. Cell alignment. For an ordered array (a) most cell attachments are either end to end 
or side to side while a disordered array (b) exhibits significant binding between poles and lateral 
surfaces. 



region. In this case, the extension, and its orientation, may be determined by the interactions 
at the boundaries which "channel" the active cells, rather than being an intrinsic collective 
property of the group of active cells. Under these experimental circumstances the boundaries 
strongly influence active cell movements. Indeed, in the physical model of Weliky et. al. [f| , 
convergent extension results only if active cells at the boundaries parallel to the elongation, 
behave differently from those at the boundaries perpendicular to the elongation. 

A subsequent and elegant experiment by Shih and Keller |It|, however, strongly sug- 
gests that, in addition, the active cells have a strong intrinsic collective mechanism driving 
their convergent extension. In these experiments a layer (essentially a monolayer) of active 
cells was excised from a frog embryo, at a stage before convergent extension had begun, 
and cultured on a uniform surface in a medium which allowed their normal development. 
Subsequently the layer showed strong convergent extension in the plane of the substrate - 
and this in the absence any plausible lateral anisotropy either in the substrate or in the 
culture medium. This behavior thus appears to be an example of "broken symmetry" so 
well known in condensed matter physics, and asks for an explanation based on collective 
behavior induced by cell-cell interactions. 

To explain this behavior by energy minimization we assume that cell-cell interactions take 
place through surface adhesion, which can be characterized by an energy per unit contact 
area. We assume that the cell rearrangements take place with negligible cell division and 
little change in cell volume, as is observed in the later stages of the above experiment. There 
is no clear understanding in the literature of the trigger for the cell elongation which initiates 
the convergent extension, and our model does not provide this. Our main assumption is 
that the adhesive energy of the contact surface between two cells will depend on how that 
surface is oriented relative to the axes of elongation of the two cells. This would be so, for 
example, were the surface density of adhesive binding sites to be different on the long side 
of a cell (parallel to the axis of elongation) from that on the short sides (perpendicular to 
the axis of elongation). We can find in the literature no compelling evidence either for or 
against this assumption. We argue here that a specific form of this assumption is a sufficient 
cause of the elongation, alignment, and intercalation resulting in convergent extension. 

We give here a two-dimensional version of our proposal since the convergent extension 
takes place in the plane of the substrate and the height of the cells does not seem to play an 
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important role. Hence we consider a collection of two dimensional cells of (nearly) the same 
fixed area. We first want to find the conditions which favor alignment. We assume a compact 
array of cells large enough that array surface effects are (for this argument) negligible. Fig. 
||]a is a cartoon of a few elongated cells in such a large ordered array of cells and Fig. 
is for a disordered array. Suppose that we can roughly distinguish, for each cell, two long 
sides (parallel to the axis of elongation) and two short sides (perpendicular). Figure §a 
shows that in the ordered array the cell-cell contact surfaces are, for the most part, either 
roughly parallel to the common axes of alignment or roughly perpendicular to that axis. We 
shall term these long-long (11) or short-short (ss) contacts since they occur, primarily, at 
contacts between a pair of long sides or a pair of short sides. The disordered array of Fig. 
0b has many contact surfaces that make intermediate angles with the now different axes 
of adjacent cells. We term these long-short (Is) contacts since they tend to occur whith 
the contact surfaces between a long side of one cell and a short side of a neighbor. If the 
energy density (per unit length) of the Is contacts is enough larger than those of and ss 
contacts then the ordered array will have the lower energy per cell (we have assumed that 
the array is large enough that we can neglect the effect of the array boundaries on the bulk 
ordering). More quantitatively, let I and s be the average long and short side lengths of 
each cell, which, for the moment, we take as fixed. Suppose that all cell-cell contacts can be 
characterized as 11, ss, or Is and that the total lengths of each type in the array are Lu, L ss , 
and L[ s . in a large array of iV cells we have 2NI = 2Lu + L/ s and 2Ns = 2L SS + L/ s (again 
neglecting array boundaries where cells do not contact other cells). Since N, I and s are 
fixed these equations provide two constraints between the three contact lengths. We assume 
that three energy densities (Ju, J ss , and Ji s ) are adequate to characterize the interactions 
at the various surfaces. Then the bulk energy of an array due to cell-cell interactions is 

E = LuJn + L SS J SS + L ls Ji s (1) 
= (2NI - L ls )J u /2 + (2Ns - L ls ) J ss /2 + L ls J ls 
= (J ls - J u /2 - J ss /2)L ls + N(U U + sJ ss ). 

This energy is an increasing function of Li s if the ordering condition: 

lis = Jis - (Ju + J ss )/2 > 0, (2) 

is satisfied. In this event ordered arrays (Li s = 0) will have lower bulk energies than 
disordered (Li s > 0) arrays. Note that condition @ is just that the Is surface tension 7^ 
be positive. 

The above argument is exact if the cells are assumed (unrealistically) to be identical 
rectangles arranged in arbitrary tesselations of the plane and is similar to that used in || in 
the cell sorting problem. For realistic cells it is a crude but plausible representation of the 
assumed anisotropy of the surface adhesion. 

The case where the adhesive energy density of a two cell contact is the product of a factor 
from each cell is interesting. For example, where a variation in the density of the binding 
sites on the cell surface causes the variation in energy. If we make the natural assumption 
that the density of adhesive bonds is proportional to the product of the density of binding 
sites on the cell surfaces in contact, then we would have in the above model Ju = —jiji, 
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J ss = —j s js, and Ji s = —jijs, where the sign is chosen make all J < when all j > 0. This 
choice satisfies the ordering condition Eq. (|l|) whenever ji and j s are positive and are not 
equal. 

In addition to Eq. ([!]) let us suppose that the // energy density is lower than the ss energy 
density. 

J u < J ss < (or ji > j s > 0). (3) 

Now the energy Eq. (|I|) of the array can be reduced by increasing the cell long side length I 
and decreasing the short side length s causing, or at least favoring, elongation of the cells. 
At equilibrium these surface effects will presumably be balanced by internal cellular forces 
opposing further elongation. 

We can also argue that Eq. (§) will produce intercalation in the direction of elongation. 
Consider the effect of the boundary on a finite array of N cells. Suppose that boundary 
cells have no adhesive energy with the culture medium. Then the expression Eq. (JJ) un- 
derestimates the array energy because it assumes all cell surfaces are in contact with other 
cell surfaces and so overestimates the contact lengths L u and L ss . From Eq. ([[]) we should 
subtract the (negative) adhesive energy that is not present at the contacts between the 
boundary cells and the surrounding medium. Figure |1| shows arrays of twelve elongated 
cells. In Fig. |l|c the array extension is at right angles to the cell elongation and in Fig. 
|l|b it is along the cell elongation. Clearly in Fig. |l]c the boundary contacts are primarily 
through short cell sides whereas in Fig. |l|b they are primarily though long cell sides. Since 
the long sides have lower (more negative) energy than the short, the energy (corrected for 
boundaries) of the configuration shown in Fig. |l|b is higher than that of Fig. |l|c. Thus if we 
start with any compact initial array of unelongated cells we expect cell motility and energy 
minimization to produce configuration of type Fig. |l]c by cell elongation, alignment, and 
intercalation parallel to the alignment. In order for these processes to produce net exten- 
sion in the direction perpendicular to alignment the effects of intercalation must dominate 
those of elongation. For a rectangular array of a large number of rectangular cells a direct 
calculation shows that this will be so, independent of the degree of elongation, and that the 
ratio of the array dimensions in the directions perpendicular and parallel to the elongation 
is Ju/ J ss - We derive these results more generally below. 

We can make the above arguments concerning surface effects somewhat more realistic 
and quantitative by the following mean field type of modeling. We assume that we have a 
large array of N elongated and aligned cells. The total energy of the array is the bulk energy 
due to cell-cell interactions plus the surface correction for the absence of cells outside the 
boundary. The bulk energy is proportional to N, or equivalently, the array area A, so we 
write it as XA, where A is the (negative) bulk energy per unit area in the aligned array. To 
model the anisotropic cell-cell interaction we assume that J depends on the angle between 
the direction of alignment, specified by the unit vector a, and the unit vector n normal to 
the contact segment between the cells (see Fig. |3]a). More explicitly, we assume that J(n-a) 
is negative, an even function (since a, —a and n, — n specify the same physical situations), 
and is minimum at na = (so that 11 interactions have the lowest energy). Figure |3|b shows 
part of a finite array of vertically aligned cells and their boundary with an external medium 
with which we assume they have no adhesive energy. To get the energy of the finite array 
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FIG. 3. Anisotropic binding. Adhesive energy at the point of contact between cells (a) is 
assumed to depend on (n • a) 2 where n is the local unit normal while a gives alignment, assumed 
common to all cells. At an interface with uniformly inert surroundings (b), missing adhesive energy 
will vary with the orientation of the surface cells, relative to the boundary. 



we must subtract from the bulk energy half the energy the boundary cells would have had 
with cells external to the array had the boundary been absent. Half, since adhesive energy 
is shared between two cells. So: 

E = XA - - j> J(n ■ a) dl, (4) 

where the integral is taken around a closed boundary. We want to minimize this over 
all closed boundaries enclosing the same area A. Alternatively we can interpret A as a 
Lagrange multiplier and find the extrema of Eq. @ over all closed curves at fixed A. To 
do this parameterize the curves as r(u) with < u < 1, and r(0) = r(l). Then, since dl = 
(x 2 + y 2 ) 1 ^ 2 du (where r = dr/du), while (n • a) = (a y x — a x y)/(x 2 + y 2 ) 1 ^ 2 and A = Jq yx du 
we can write the energy as Jq 1 £(r, r) du with £(r, r) = Xyx — J(n • k)(x 2 + y 2 ) 1 ^ 2 /2. The 
extremal curves are solutions of the usual Euler-Lagrange equations for C and are degenerate 
with respect to translations in the x-y plane. This degeneracy gives rise to two first integrals 
and two constants of integration (which we choose to be zero), which fix the position of the 
extremal curve. The integrated equations have the form 

2Ar = aJ'(n ■ a) + n[J(n • a) - (n • a) J'(n ■ a)], (5) 

where J' is the derivative of J. The dot product of Eq. [5] with f gives the Wulff condition 



TTJ. Equation (|Sp is two coupled first order differential equations whose solutions depend 
on the particular choice of the function J. We have not been able to find complete analytic 
solutions for any interesting choice of J but some properties of the solutions can be found. 
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First we note that for a = 0, or equivalently J = constant, the solution is a circle of radius 
J/(2A). Secondly, the turning points of any solution curve are where d(r ■ r)/du = 0. Now: 

A^-{r t) = 2A(r-r) = (r • a) J'(n • a), (6) 

where we have used Eq. ([5]) and that r • n = for any curve. From Eq. ((^) we find two 
types of turning points. 1) At r • a = 0, that is, where the boundary is perpendicular to the 
alignment so that n = ±a. For any simple closed curve this condition will be satisfied at 
two points on the curve. Because J is even and J' is odd we have from (||) that these two 
points lie at ±J(1)/(2A) on the line through the origin and parallel to a. 2) At n ■ a = 
where J' = and the boundary is parallel to the alignment. For these ([5]) shows that 
2Ar = nJ(0), thus there are turning points at ±J(0)/(2A) along a line through the origin 
and perpendicular to a. If we let D± and Du be the distances between the turning points 
aligned respectively perpendicular and parallel to a then the aspect ratio of the boundary is 

DJD\\ = J(0)/J(1). (7) 

If J(l) < J(0) < then D±/D\\ > 1 and the elongation is in the direction perpendicular 
to the alignment, as is observed in convergent extension. From Fig. [5|a we see that J(0) 
corresponds to our previous Ju while J(l) corresponds to J ss . 

We have also studied the minimization of the energy functional Eq. (|^) numerically for 
the case where J is chosen to be a gaussian function. We approximate the boundary curve 
by a polygon of at least 100 sides and use an iterative process that moves down the energy 
gradient at constant area. We have started from many initial configurations, all of which are 
simple closed polygons. The final boundary curve is always the same and with the correct 
aspect ratio Eq. (|?p. This iterative process could also be viewed as a model for the dynamics 
of convergent extension. Indeed, with the addition of additive random forces, the method 
would be essentially a Langevin dynamics for the evolution of a highly dissipative system. 

In conclusion, we can understand convergent extension as an energy minimization pro- 
cess, provided the cell-cell adhesive energy has a certain kind of anisotropy. This single 
simple property is sufficient cause for the cell extension, alignment, and intercalation in 
the direction of alignment, that are the characteristics of convergent extension. We have 
characterized the anisotropy required [Eq. (0) and Eq. @]. We have shown that the final 
aspect ratio is independent of the initial configuration and have shown how it depends on 
the anisotropy Eq. (|7|). 

We believe our arguments are plausible but realize that they are not conclusive. Our 
modeling neglects many degrees of freedom associated with cell shape and arrangement, 
which we think, but cannot prove, are not crucial. Our procedure of separately minimizing 
the bulk and surface energies is accurate only for a large array of cells. We do not see much 
possibility of doing a lot better by purely analytic methods. We have initiated simulations of 
convergent extension, using the Potts model and Metropolis dynamics methods of references 
H and [0, with anisotropic adhesive energies of the type described in this paper. The use 
of anisotropic adhesive energies introduces technical difficulties in that the energy becomes 
non-local on the scale of the size of a cell, which considerably increases the simulation time. 
Nevertheless we believe the simulations will eventually substantiate our conclusions. Even 
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so, the more difficult question of whether this explanation of convergent extension is correct 
remains. Experiments that probe the possible anisotropy of cell adhesive energy would be 
useful, as would experiments that show the final configuration is largely independent of the 
initial configuration. 
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